Topological r/— SQUID with a frustrating tt— junction for probing the Majorana 

Bound State 
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Majorana Bound States are predicted to appear as boundary states of the Kitaev model. Here 
we show that a tt— Josephson Junction, inserted in a topologically non trivial model ring, sustains 
a Majorana Bound State, which is robust with respect to local and non local perturbations. The 
realistic structure could be based on a High Tc Superconductor tricrystal structure, similar to the 
one used to spot the d-wave order parameter. The presence of the Majorana Bound State changes 
the ground state of the topologically non trivial ring in a measurable way, with respect to that of a 
conventional one. 



After the pioneering theoretical proposal by Kitaevjlj, 
Majorana Fermions (MPs) have been predicted in a wide 
class of low-dimensional solid state devices. Being neu- 
tral excitations in Fractional Quantum Hall systems or 
hybrid superconducting devices, MFs are highly attrac- 
tive for quantum computing gates, as well as for fund- 
mental reasons. Despite the considerable theoretical and 
experimental efforts ,2'], challenges still remain before a 
real solid state device can be realized, allowing for isola- 
tion and manipulation of MFs. Among the promising en- 
gineered systems, there are quasi one-dimensional super- 
conductors in contact with topological insulators (TIs) [s^ 
or quasi one-dimensional systems with strong spin orbit 
interactions 043' helical magnets [1] and other materi- 
als [gl-fisl . However several issues have to be convincingly 



solved [mSI and the recently announced experimental 
realizations 17 - 1^ still arise excitement and debate (2^ . 

Herewith we provide evidence that the richer function- 
alities of High Tc Superconductors (HTS), hybridized 
with topologically non trivial materials, can be used to 
monitor unambiguously the emergence of MFs in the 
ground state (GS) of a SQUID Josephson ring. 

Epitaxial growths of high Tc cuprates, matching three 
differently orientated crystals, has provided in the past 
the evidence for the d — wave symmetry of the su per con- 
ducting order parameter (OP) in HTS materials [2l|. In 
Ref. 2l| the experiment has been properly designed so 
that the ring contains an odd number of 7r-junctions and 
has a total critical current (/c) and an inductance (L) to 
guarantee IcL » 4>o (here (po = hc/2e is the supercon- 
ductive flux quantum). Its GS spontaneously breaks the 
time reversal symmetry with a current flowing in the ring 
and generating a spontaneous fractional flux </> = ±0o/2 
that can be measured by scanning SQUID microscopy 

The Kitaev model offers a very rich playground to 
study topological non trivial states. Inspired by the con- 
cept that any superconducting loop with an odd num- 
ber of TT junctions has a frustrated GS0 (we call it 
"frustrated tt— ring" (FttR ) in the following), we study 



a closed loop " Kitaev chain" [ij , with a weak link having 
a jump of = TT in the phase of the OP (see FigH] for 
a sketch). Real life realizations of such a system could 
be a semiconducting wire with strong spin-orbit coupling 
(e.g. InAs or InSb), with spin polarizations splitted by a 
Zeeman magnetic field, or, alternatively, the edge state 
of a quasi two dimensional TI flake, such as Bi2Se3 or 
Bi2Te3. However, in this letter, we leave aside the de- 
tails of the engineering of such devices and we focus on 
the model only. In particular we explicitly calculate the 
various branches in the energy dispersion of the Andreev 
Bound states En{ip) of the model ringl'], finding a cross- 
ing of states of opposite fermionic parities V at zero 
energy, which signals unambiguously the presence of the 
Majorana Bound State (MBS) in the spectrum. 

As shown in Ref. it junctions obtained by prox- 
imity induced superconductivity from nodal d-wave su- 
perconductors, can display a MBS at a phase difference 

= 0, rather than ip = n. Here we show that a cor- 
responding state exists in the ring geometry depicted in 
Fig. [T] and is robust with respect to residual interactions, 
paving the way for an unambiguous measurement of its 
presence. 

In general, the free energy describing a Josephson rf- 
SQUID ring shows minima at a phase difference ip which 
are directly related to a measurable spontaneous flux. In 
the conventional FttR geometry, there is quite a high bar- 
rier separating the two minima at ±(/)o/2, which freezes 
the frustrated system in a time reversal symmetry bro- 
ken GS. In our model, we find two minima close to ±0o/2 
as well, as in the conventional FttR . However, they cor- 
respond to different parities. While the isolated model 
Hamiltonian conserves V, this symmetry is rather un- 
likely to be mantained in the real life. In fact, impuri- 
ties, trapping and releasing charges in the background, 
cannot be avoided. These events would correlate to flux 
quanta entering or leaving the ring, thus inducing tun- 
neling of the system between the two minima, so that the 
expected average flux associated to the GS is ($) = 0. 
Measuring zero flux at a FttR in presence of topologically 
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FIG. 1: (color online) Sketch of the Kitaev chain with a 
TT— junction at the top point U. Red (blue) circles represent 
7a(s) type MFs. Yellow ellipses signal strong coupling. 



protected states would be the smoking gun evidence for 
the existence of the MBS in this device. It is a weird case 
that the breaking of the discrete symmetry (V) enforces 
the time reversal symmetry to be restored. 

Model Hamiltonian 



We start considering a N-sites Kitaev chain of unitary 
lattice spacing and full length L = 2N, with real inter 
site hopping t. When folded in the shape of a ring the 
system displays mirror symmetry across the vertical line 
connecting points U and D between the top and the bot- 
tom (see Fig. [1]). 

Our device can be described as two N-sites Majorana 
wires (left {£) and right (r)), coupled at the top of the 
ring U by the weak electron tunneling of energy F. At 
chemical potential = 0, in the presence of an electro- 
magnetic vector potential A, the gauge invariant Hamil- 
tonian reads as H — Hi -\- H^, with: 
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Here a labels the £ and r side of the ring, Caj are spinless 
Dirac fermions at the site j and A is the superconducting 
pairing, generating an effective p-wave superconductivity. 
The phases gaj acquired in the hopping between the j site 
and its nearest neighbor and the gauge invariant phase 
(f>aj are defined as: 

g paj + l ^ 

y^^ A-di, (2) 

2e - ^ 
0aj ^0^j~— A-dl, (3) 
tic Jii 

where 9aj is the phase of the superconducting order pa- 
rameter. 

At the top point U of the ring (see Fig.l) there is a 
tunnel junction which allows for the fractional Josephson 
coupling: 

Hu = -F (^4^CrN + /i.e.) , (4) 

(F << i). For sake of further investigations, we explicitly 
consider also the hopping term at the bottom point in the 
ring D, where the £ and r chains are matched: 

Hd = -ut [c\^Cri + h.c] . (5) 



Here we will keep the dimensionless parameter u (which 
may be complex) as a variable, to discuss also the limiting 
case of w = 0, which corresponds to the ring cut at D, 
with open ends. 

The spinless Dirac fermions can be expressed in terms 
of two species of Majorana fermions at each site of the 
ring lA/Bj-' ^^^^ th^t- 

^g^.=c„,e^^-/^+4■e-'^-/^ (6) 
73,--*(ca,e^'-/2-4.e-^-/2) . (7) 

A TT— Josephson Junction requires A having opposite 
signs at U, between £1 and rN . In the gauge in which 
A is real, the OP A has to vanish somewhere along the 
ring and we choose this point to be D with no loss of 
generality. As a first step, to make the approach as sim- 
plest as possible, deep in the topological phase, we will 
adopt the Kitaev approximation, |A| =t [jjall along the 
chain and we choose A = i in the £ region and A = — t 
in the r region of the ring. Thus, the chain Hamiltonian 
becomes: 

He + Hr = 2 E bB,7Aj+i - 7aj7b,+i] ■ (8) 
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Pictorially, this kind of hybridization can be represented 
as in FigUl Blue (red) circles represent B (A)-type MFs. 
The yellow ellipses denote effective strong coupling be- 
tween nearest neighbor MFs. Were Eq.Q the full Hamil- 
tonian, the t and r chains would dimerize with opposite 
phases. In addition, four unpaired MFs appear in Fig.l: 
two (the red/ A ones) located at U and two (the blu/B 



ones) at D. 

To account for the extra interactions F and ut, fol- 
lowing Ref. [2^ we refermionize the Hamiltonian by in- 
cluding Hu -f Hd and by rearranging the MFs at the 
boundaries. Three effective Dirac Fermions, c?Aj i^b, dend 
are required, located at the U weak link, and three more 
ones, Ja^ fs, fend, at the D boundary, according to[26|: 
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where 
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is the phase difference at the U weak link. The to- 
tal energy only depends on the flux threading the ring, 
$ = § A ■ \n units of 0o- This Bogolubov-de-Gennes 
(BdG) Hamiltonian i?e// changes sign under the opera- 
tion S of exchanging particles with holes and conserves 
the Fermionic parity V = (^—\)'i^ifiSi+didi_ can be 
shown that the low energy spectrum only depends on 
the first term in the second line, which involves the Dirac 
fermion d^nd = {iai + * 7a7v) /^i and on terms involving 

fend, find "^1*^ fend = (tbi + « TsAf) /2, which cau be 
obtained by perturbation theory from the first term of 
the first line. In addition, residual interactions, not in- 
cluded in Heff, can account for finite size effects. In 
particular an extra coupling z^F cx {a = i, r) arises 
from realistic long range interactions between the edge 
MFs at each chain, when the simple t = |A|[l[ limitation 
is relaxed. 

We are led to the minimal 4X4 Hamiltonian in the 
Majorana representation, just involving the four relevant 
MFs: 
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(11) 



in the basis [7^^, 7^^, 7^^, Ybi]^ with u' = ut/T. 

Ring cut at the bottom point : u — Q. 
When the control parameter u is set to zero, the ring is 
cut at the bottom point D and the system is equivalent 
to a linear topological tt junction with phase difference Lp 
and open ends. In this case, if the finite size couplings Zq, 
are neglected (see Fig. [5] left panel), there is a crossing at 
zero energy and zero flux due to the MFs at the U point, 
signalling a change of parity in the GS when the flux 
changes sign. Together with it, two dispersionless zero 



energy modes appear, corresponding to the dangling MFs 
at the open ends, However, the system is expected to be 
unstable with respect to finite size interactions described 
by the couplings. Indeed as soon as one turns z^ on, a 
gap opens and the zero energy MBS disappears (see Fig. 
Elb). 
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FIG. 2: (color online) Andreev Bound States as a function of 
the magnetic flux for an open ring {u' = 0). Left panel) Finite 
size corrections not included zj = Zr = 0. Right panel) When 
finite size corrections are included zi = Zr = 0.2 there is a 
splitting of the zero energy Andreev Bound State. 



TT—ring configuration: u ^ 0. 
In the ring configuration (u 7^ 0), the situation is quite 
different (see Fig. [3]). The zero energy MBS is always 
present, no matter how strong the finite size effects 
arc. In the non-symmetric case (zg ^ z^), the location 
of the crossing occurs at non zero flux. With increasing 
of u' , the flux of the crossing point drifts towards ip — 
(see Fig. [3] bottom panels). For the physically relevant 
case u' >> Zg, Zr, just the crossing Andreev bound states 
survive at low energy and the dispersion turns out to 
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be approximately symmetric. The dispersion of the two 
crossing low lying energies tends to 

E±^T sin (I) (2ne„d - l)p = ±r sin (|) . (12) 

With the labeling of the average {■■■)v we denote the 
fermionic parity of the GS. The minimum is for tp « ±7r 
(i.e. flux ±(/io /2), depending on the occupancy observable 
nend = d\^^dend of the MBS located at the Josephson 
Junction. 



u'=0.i Zj = b.2 zj = 0.3 
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FIG. 3: (color online) Andreev Bound States as a function 
of the flux for the n ring. Left panels) No asymmetry z — 
z( = Zr. Right panels) Asymmetrical case: ze ^ Zr- An 
approximately symmetric spectrum is recovered for sizable 
u' , independently of the values of the asymmetry parameters 
Za- In the closed ring geometry the zero energy Majorana 
Bound State is always present. 



Model free energy and stationary conditions of the 
rf -SQUID. 

We have shown that, in the topologically non trivial 
TT— ring structure, an unpaired zero energy MBS exists 
and it is robust with respect to pertubations. This shows 
up as a crossing of the particle and hole excitation dis- 
persions at flux close to zero. Our tt— ring modelizes a 
topologically non trivial r/— SQUID device and we now 
argue that the MBS characterizes in a measurable way 
the stationary conditions of the device. 

If the ring has a small diameter, so that its self- 
inductance L cannot be disregarded, the free energy is 
a function of the circulating current and of an external 
flux (jjext which may be intentionally added. Its simplest 
form, arising from EafT^is: 



Fp(/,, 



1 



-LI' 



<1>qIc 



2 27r \<f>o 

(13) 

The free energy F-p (/, ipext — 0) at zero external flux is 
plotted in Fig. 2^). The first and second minimum, be- 
longing to the same V , differ in phase by ~ 2tt. Fig|3)D 



reports the change in shape of the free energy for one 
single V, at different applied fluxes ipext- 

The similarity with the conventional YBCO FttR is 
only superficial. At first sight Fig. |4)d could report the 
free energy plot of a conventional FttR . By fine-tuning 
the external flux the two energy minima can be made de- 
generate. If Vis strictly conserved, this would be the end 
of the story and no difference would arise. However, in 
the topologically non trivial case, there is a correspond- 
ing set of curves belonging to the other V. A switching 
between the two different minima e.g. at (jiext = (see 
Fig. 2^) not only requires a fiux quantum entering or 
exiting the ring but a simultaneous change of V . 
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FIG. 4: (color online) a): free energy vs LI/(j)o at zero ex- 
ternal magnetic field for (f)ext = {u — and ze = Zr =) ■ 
Full (dashed) line belong to the two different fermion parity 
V . h): Corresponding free energy for different values of the 
external magnetic flux and one single parity. The curves are 
displaced in energy for clarity 



A tool to fix the parity could be a side quantum point 
contact (QPC) controlling the charge tunneling. Addi- 
tion of an electron on the ring would suddenly require the 
switching of the whole device between the two possible 
GSs, corresponding to a jump in the trapped sponta- 
neous flux. By contrast, a conventional FttR is expected 
to be widely insensitive to the in-out tunneling of in- 
duced charges. Operating with a side gate on the QPC 
allows to distinguish a topological non trivial ring from 
a conventional one. 

However, in real life, the job of fixing T^appears to be 
rather hard. The ring is not expected to be isolated: 
background impurities could provide charge noise, by re- 
leasing or capturing charges. For a system open to the 
environment, the energy spectrum of the isolated ring 
''looses meaning and a description of the state of affairs 
in terms of the statistical density matrix p{t) is required. 
The latter accounts for the transitions, with absorption 
and emission of the energy between the two parities and 
with simultaneous switching of the flux. Under these 
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conditions, the expectation value of the observable cor- 
responding to the flux, is likely to average to zero: 

lim ($) = Tr {p $} ~ (14) 

t— >cx> 

Hence, we conclude that if we prepare the system in one 
of the two wells with initial flux ^ ±(/)o/2 and we leave 
the system to evolve freely by turning off the external 
flux, a measurement of the flux threading the ring after 
relaxation has taken place gives a null result. This is 
in striking contrast with a conventional tricrystal YBCO 
structure, which would not abandon the state support- 
ing the spontaneoiis half flux quantum, and stays as a 
distinctive fingerprint of the presence of MFs. 
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